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Background : 

The subject grant no. NAG 3-506 was initiated in December 1983 
and ended on January 1990. During this period, with the support 
of this grant several small research projects were carried out 
and a number of undergraduate and graduate students of the 
Department of Materials Science and Engineering, University of 
Cincinnati were trained at NASA Lewis Research Center under the 
supervision of various scientists. The research work was 
conducted at the University of Cincinnati by the some of the 
students under the supervision of the P.I. The following 
students were involved in these activities: 

1. Walter W. Milligan - Senior Project thesis, 1983-84 - (NASA 
Advisor : Dr. Robert C. Bill) 

2. Kevin Hemker - Senior Project thesis, 1984-85 (NASA Advisor: 
Dr. Robert V. Miner) 

3. Dan Thoma - Senior Project thesis, 1985-86 (NASA Advisor: Mr. 
Thomas Glasgow) 

4. Pete Kantzos - Summer student at NASA-Lewis Research Center, 

1985 (NASA Advisor: Mr. Jack Telesman and Dr. Marv Hirschberg) 

5. David Lee - Summer student at NASA-Lewis Research Center, 1985 
(NASA Advisor: Dr. Hugh Gray) 

6. Frank Ritzert - Summer student at NASA-Lewis Research Center, 

1986 (NASA Advisor: Dr. Hugh Gray) 

7. Vi jay Ramakrishnan - Graduate student training at NASA-Lewis 
Research Center, Winter 1989 (NASA Advisor: Dr. Brad Lerch) 

8. Thomas E. Fetsko - Graduate student training at NASA-Lewis 
Research Center, 1988-90 (M.S. degree awarded: 1990; NASA 
Advisor: Dr. Ram T. Bhatt). 

9. Monika Ditmars - Graduate student at the University of 
Cincinnati, 1984-86 (M.S. degree awarded: 1987) 

10. Jun Sonu - Graduate student at the University of Cincinnati, 
1984-88 (Ph.D. degree awarded: 1988) 

Results from these research work have been presented at various 
NASA meetings, National and international conferences and some of 
these work have also been published as journal articles. This 
report summarizes some of the important results from two of the 
very recent projects completed under this program: 


Multiaxial fatigue damage evolution in waspaloy: 

ABSTRACT 

Strain controlled torsional and biaxial (tension-torsion) 
low cycle fatigue behaviour of Waspaloy was studied at room 
temperature as a function of heat treatment. Biaxial tests were 
conducted under proportional (when the axial and torsional strain 
cycles are in-phase) and non-proportional (when the axial and 
torsional strain cycles are 90° out-of-phase) cyclic conditions. 
The deformation behaviour under these different cyclic conditions 
were evaluated by slip trace analysis. For this, a Schmidt-type 
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factor was defined for multiaxial loading conditions and it was 
shown that when the slip deformation is predominant, non- 
proportional cycles are more damaging than proportional or pure 
axial or torsional cycles. This was attributed to the fact that 
under non-proportional cyclic conditions, deformation was through 
multiple slip as opposed single slip for other loading 
conditions, which gave rise to increased hardening. The total 
life for a given test condition was found to be independent of 
heat treatment. This was interpreted as being due to the 
differences in the cycles to initiation and propagation of 
cracks . 


1. INTRODUCTION 

Aircraft engine components are often subjected to complex 
multi-axial cyclic stresses. Plastic deformation behaviour under 
multi-axial stress state is different from that observed under 
uniaxial conditions (1-4) . In this paper the plastic deformation 
aspects of multiaxial fatigue is considered. Waspaloy was chosen 
for this study since its uniaxial fatigue behaviour was 
extensively studied in an earlier study (5) . 

Lerch, Jayaraman and Antolovich (5) examined the uniaxial 
low cycle fatigue (LCF) behaviour of Waspaloy. Two different 
heat treatments and therefore, two different microstructures were 
examined, namely heat treatment A to produce a microstructure 
consisting of fine grain size (ASTM G.S.#9) with large gamma 
prime prcipitate particles (about 900 A diameter) and heat 
treatment F to produce a microstructure consisting of coarse 
grain size (ASTM G.S.#3) with small gamma prime precitate 
particles (about 50-80 A diameter) . Although under strain 
controlled fatigue testing both these microstructures gave 
relatively similar total lives, for heat treatment A most of the 
life was spent in crack initiation and the plastic deformation 
was uniformly distributed by the Orowan looping mechanism. As 
opposed to this, for heat treatment F, most of the fatigue life 
was spent in propagation of cracks initiated earlier, and the 
deformation was localized in slip bands formed due to precipitate 
shearing mechanism. 

In this paper, again Waspaloy with the above two 
microstructures was tested for the torsional and biaxial 
(tension-torsion) LCF behaviour. Biaxial fatigue studies where 
the tension and torsional cycles were in-phase and where the two 
were 90° out-of-phase were conducted. The main purpose of these 
studies was to illustrate the slip and deformation processes, and 
fatigue fracture in the different types of biaxial fatigue 
conditions. 

2. MECHANICS OF BIAXIAL FATIGUE 

2.1. Schmidt factor calculations for slip 

Biaxial LCF testing in this study was accomplished by 
cyclically applying torsional and axial strains to a tubular 
sample. Both cycles were applied at the same frequency. A 
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biaxial test is said to be in-phase, or proportional, when 
a /r or e /y = constant. It is said to be out-of-phase, or non- 
proportional when there is a phase shift between the axial and 
torsional cycles. Figure 1 demonstrates the difference between 
in- and 90° out-of -phase fatigue cycles. Note that purely axial 
and purely torsional tests are proportional. 

When a sufficient stress or a combination of stresses is 
applied to a metallic material, plastic deformation occurs. On a 
microstructural level, this deformation is caused by the motion 
of dislocations (slip) . In face centered cubic materials, slip 
tends to occur usually in the twelve possible octahedral slip 
systems (denoted by (111}<110>) . For any given. state of stress, 
slip is most likely to occur on those slip systems with the 
greatest resoved shear stress. 

In the most general case the resolved shear stress can be 
determined as follows (6) : 


a XX 

7 xy 

7 xz 
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where <f> is the unit vector in the direction of the slip plane 
normal, 5 is the unit vector in the slip direction and a ^ j is a 
stress tensor describing the state of stress at the point under 
consideration. For the case of the simple uniaxial stress, this 
equation reduces to a more familiar form as follows: 

t j. o Cos ^ Cos 6 . . * ( 2 ) 

which is the Schmidt equation. The term Cos <f> Cos 6 is the 
Schmidt factor and has a maximum possible value of 0.5. 

For the more complex biaxial case also, Schmidt-type factors 
can be determined. For tension-torsion loading the Schmidt 
factor has been shown to be (7) : 


7 RSS 


{a 2 + 4 t 2 ) 


... (3) 


where a = applied axial stress, and 
t = applied shear stress. 

The biaxial Schmidt factor also has a maximum possible value 
of 0.5. 
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2.2. Planes of Maximum Shear 

For any given stress state, the axes of principle stresses 
and the planes of maximum shear stress can be uniquely determined 
by using standard tensor algebra (6) . It can easily be 
demonstrated (7) that for a cylindrical or tubular body stressed 
by any combination of axial and torsional loads that two of the 
axes of principle stresses lie on (or parallel to) the plane 
tangent to the surface of the body at any point on its 
circumference. The third principle stress is normal to the 
surface and always has a value of zero. The maximum shear stress 
acts on a plane whose normal bisects the angle. between the 
directions of the largest and smallest principle stresses. The 
planes of maximum shear is therefore perpendicular to the surface 
of the body at all points (figure 2) and their orientation with 
reference to the sample axis has been shown to be (8) : 

tan 29 = — • • • (4) 

2r 


where 0 is the angle between the specimen axis and the plane of 
maximum shear. 

Under pure axial or pure torsional cyclic loading the 
orientations of both the principle axes and the maximum shear 
planes remain fixed with respect to the sample during the course 
of each cycle. For pure axial loading the maximum shear planes 
can be shown from equation 4 to be at angles of 45° and 135° to 
the specimen axis and for pure torsional loading these planes are 
parallel and perpendicular to the sample axis. This means that 
the Schmidt factors for various slip systems do not vary during 
fatigue testing unless the crystal itself rotates. In other 
words, the resolved shear stress will always be greatest on one 
slip system (or possibly several systems if the Schmidt factors 
on these are identical) . 

Figures 3a and b demonstrate how the resolved shear stress 
on the twelve possible slip systems varies during on axial and 
one torsional fatigue cycle, respectively. The applied load was 
varied sinusoidally for these examples and the cycle begins and 
ends at zero load. The crystal in these examples was arbitrarily 
chosen to have a single slip orientation of [158] for the 

specimen axis and [7"31] with respect to the surface normal. 

Notice in these figures that one slip system always has the 
maximum shear stress at all times during the cycle. 

During a biaxial fatigue cycle, the principle axes and 
therefore maximum shear planes maintain a fixed orientation if 
and only if the loading is proportional, i.e., o/t = constant at 
all points during the fatigue cycle. For example, for the test 
conditions of this study, the orientation of the maximum shear 
planes were determined using equatin 4 to be at fixed angles of 
30° and 120° with reference to the sample axis. Figure 3c shows 
how the resoved shear stress varies on the twelve slip systems 
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for a proportional biaxial fatigue cycle. Notice again, that 
only one slip system has the largest resolved shear stress value 
at all times during the cycle. 

If the loading is non-proportional, the principle axes and 
therefore the maximum shear planes will rotate with respect to 
the crystal by 180° in the course of each cycle. Figure 3d shows 
how the resolved shear stress on the twelve slip systems varies 
during a 90° out— of-phase cycle. Notice that now several slip 
systems share the privilege of having the greatest resolved shear 
stress during various portions of the cycle. Under these 
circumstances, extensive multiple and/or cross slip is highly 
possible. 

In this paper, it is shown that the above analysis is very 
consistent with the observed slip and hardening behavior of 
Waspaloy subjected to torsional and biaxial (both proportional 
and non-proportional) fatigue. 

3. EXPERIMENTAL 

3.1 Materials and Heat Treatments 

A Waspaloy billet was supplied by Cyclops Cytemp Specialty 
Steel Division. The chemical composition of the billet is listed 
in Table I. The 55cm X 50cm X 5cm billet was hot rolled at 
1230°C to approximately 100cm X 50cm X 3.2cm. The final hot 
rolling was done in the 100 cm direction. The plate used in the 
previous investigation (5) was cross rolled to its final 
dimensions. As will be shown later, this final hot rolling did 
not introduce any significant anisotropy in the grain structure 
of the billet. 

Specimens were heat treated in air, in a resistance-type 
furnace. The heat treatments were the same as those used in the 
previous investigation ( 5) . Details about the heat treatment and 
the microstructures can be found in reference (5) . For the sake 
of completeness these are summarized in Table II and discussed 
briefly in Section 4. 

3.2. Mechanical Testing 

3.2.1. Specimen Configuration 

Longitudinal LCF specimens with a hollow cylinder 
configuration were machined from the billet with their long axis 
parallel to the final rolling direction. An outer diameter (o.d.) 
to wall thickness ratio of 10:1 was chosen to facilitate stress 
calculations. With this ratio, the specimen can be considered to 
be thin walled and the stresses can be assumed to remain the same 
through the wall thickness. Figure 4 shows the specimen 
geometry. The central bore was gun-drilled and honed to a 16 
micron rms (root mean square) finish. The outside surface was 
low stress ground and polished to a 8 micron rms. In addition, 
all the specimens were also electropolished on the outside 
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surface and examined under a low power optical microscope for any 
surface scratches or damage. 

3.2.2. Fatigue Testing 

Fatigue tests were performed at room temperature on a closed 
loop servohydraulic biaxial (tension-torsion) MTS machine 
interfaced with an LSI 11-23 computer. This computer was used 
for test control, data acquisition and analysis. All tests were 
strain controlled and fully reversed. Strains were detected by a 
previously calibrated MTS biaxial extensometer . The calibration 
procedures are detailed in reference (7). The extensometer is 
capable of measuring axial displacements and torsional rotations 
simultaneously. The load cell is capable of measuring both axial 
loads and torques. Therefore the axial stress and strains could 
be determined from the axial load and displacement. Since the 
torsional strain varies linearly with radius, whether the 
deformation is elastic or plastic, c = rv/L, where c is the 
torsional strain, v is the angle of twist, r is the radius and L 
is the gage length of the sample. Torsional stress is however, 
somewhat more complex to determine, because the relationship 
between stress and radius can not be defined in the plastic 
region. For this purpose, a thin-walled tube with an o.d. to 
wall thickness ratio of atleast 10:1 was chosen for the specimen 
configuration. Then assuming that the torsional stress is the 
same through the thickness of the tube (8) 


3T 

r = 

2 * ( r o 3 " r i 3 ) 


... (5) 


where T is torque and r Q and r^ are the outer and inner radii of 
the tube respectively. 

Torsional, in-phase tension-torsion and 90° out-of-phase 
tension-torsion tests were conducted for both heat treatments. 
Since plastic strain controlled (constant De_) uniaxial LCF tests 
were conducted in the earlier investigation Y 5 ) / about half of 
the tests in the present work were done under constant plastic 
strain range. In order to accomplish this, the total strain 
range was continuously adjusted during the strain controlled 
tests. This technique worked well for the torsional and the in- 
phase biaxial tests. However, in the 90° out-of phase biaxial 
tests, the definitions of axial and torsional plastic strains are 
not clear. For this reason ten of the twenty three tests were 
run under total strain control. As shown later, the test results 
from these two different types of tests showed similar trends and 
therefore could be analyzed by the same analytical methods. All 
fatigue tests were run to failure, which was defined as a sharp 
drop in the tensile load amplitude (for biaxial tests) or the 
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shear load amplitude (for torsion tests). Physical observation 
showed that this sharp drop corresponded with the development of 
a rapidly propagating macrocrack. No special effort was made to 
detect the moment of crack initiation. 

3.3. Me t a 1 1 ogr aphy 

3.3.1. Optical Metallography 

The following three different optical metallography 
techniques were used to examine selected test specimens: 

(i) Macrophotography of the crack, 

(ii) Microphotography of the gage surface, without any surface 
preperation of the test specimen, and 

(iii) Microphotography of the mounted, polished and etched 
sections of the tested samples. Details of the procedures can be 
found in reference (7) . 

3.3.2. Scanning Electron Microscopy (SEM) 

Fracture surfaces of selected specimens were examined by SEM 
to characterize the fracture mode and propagation behaviour. The 
gage surface was also examined for secondary cracking and slip 
band formation. 

4. RESULTS AND DISCUSSION 

4.1. Material Characterization 

4.1.1. Grain Size 

The grain sizes were determined by the Heyn intercept method 
given in ASTM Standard E112 (10). Heat treatment A (with a 
solutionizing temperature of 1010°C) was found to have a grain 
size of ASTM size 7 (approximately 32 microns) . No grain growth 
took place during the 1010°C solution heat treatment. Heat 
treatment F (with a solutionizing temperature of 1030°C)had a 
grain size of ASTM size 2 (approximately 180 microns) . The grain 
sizes of the two heat treatments used in this investigation are 
slightly larger than those used in the earlier work (5) and this 
is due to the difference in the grain sizes of the as-received 
billets in the two studies. In reference (5) , the starting billet 
had an ASTM size of 9-10 as opposed to an ASTM size of 7 in this 
work. As shown later, this small difference in the grain sizes 
did not significantly affect the results. Thus, as in the 
previous study, the two heat treatments were designed to give a 
fine and a coarse grain structure. 

4.1.2. Gamma Prime Phase 

The aging heat treatments (for both A and F) were designed 
to produce large and small gamma prime precipitate particles 
respectively, which in turn can affect the deformation mode 
(i.e., either looping or shearing) and hence the damage 
mechanisms. Accordingly, heat treatment A produced a duplex 
gamma prime structure with particle sizes of 900 and 3000 A. The 
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larger particles were unsolutionized and therefore grew during 
aging. Heat treatment F produced a small gamma prime structure 
with particles in the size range of 50-100 A. High resolution 
darkfield transmission electron microscopy (TEM) images of the 
precipitates for both heat treatments are shown in figure 5. 

These results are in agreement with those previously reported (5) . 

4.1.3. Carbides 

As in the previous study (5) , the fine grained material 
(heat treatment A) showed discrete carbide particles at or near 
grain boundaries while the coarse grained material (heat 
treatment F) showed an absence of grain boundary carbides. Figure 
6 shows SEM photographs of typical grain boundaries from these 
heat treatments . 

4.2. Low Cycle Fatigue Test Results 

All the fatigue test results are summarized in Table HI* 

The experimental parameters such as A«, A 7 , etc., were determined 
at N f /2, where N f is the number of cycles to failure which, as 
indicated earlier, was determined by the sudden drop in the axial 
load (for biaxial tests) and torsional loads (for torsion fatigue 
tests). The analyses of life from these test results are 
discussed in another paper (11) . Here, only the deformation 
characteristics, the cyclic hardening behaviour and the fracture 
characteristics will be discussed. 

4.2.1. Hysterisis Loops 

Figures 7a-e show typical hysterisis loops (stress vs. 
strain plots for one cycle) from torsional, biaxial in-phase and 
biaxial 90° out-of-phase tests. The biaxial tests have two 
hysterisis loops, one for axial and another for torsional. The 
torsional test loops and the biaxial in-phase test loops are very 
similar to the ones obtained in uniaxial tests. These loops 
characteristically show a linear elastic region followed by a 
non-linear plastic region, indicating that there are periods of 
time when the material experiences only elasic strains (or 
stresses) . As opposed to this, the biaxial 90° out-of-phase test 
loops (figures 7 d and e) bear little resemblance to their in- 
phase counterparts. Similar loops have been reported by other 
researchers (1) and predicted by constitutive models (12) . The 
unusual shape of the loops can be attributed to the fact that at 
no time does the material undergo elastic deformation. Notice 
that there are still linear regions in these loops. These can 
not however be considered elastic regions, since their slopes do 
not equal the elastic modulus (or the shear modulus) of the 
material. For both the torsional and biaxial in-phase loops, one 
can easily identify the plastic strain ranges by conventional 
methods from the loop widths. However, for the 90° out-of-phase 
tests, because there are no elastic regions, the loop width is 
only remotely related to the plastic strain. For this reason, 
the plastic strain control method was abondoned in favour of a 
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total strain control method for the remaining ten tests. These 
are appropriately indicated in Table III. 

4.2.3. Cyclic Hardening Behavior 

Figures 8a— e and figures 9a— e show the cumulative glide 
paths (stress vs. number of cycles) for both heat treatments A 
and F respectively. Experimental results from pure torsion and 
biaxial, proportional and non-proportional fatigue tests are 
presented in these plots. Since both axial and torsional 
stresses were monitored during biaxial tests, the corresponding 
cumulative glide plots are shown. As seen in these figures, the 
general hardening-softening behaviour is similar to that seen in 
the uniaxial LCF studies (5), namely, heat treatment F showed 
fairly rapid hardening and softening to failure. In contrast, 
the heat treatment A showed slow hardening and little softening. 
Non-proportional biaxial fatigue cycles tended to harden more 
than proportional biaxial fatigue cycles or torsional fatigue 
cycles. It will be shown that this additional hardening is the 
result of multiple slip. The relative hardening behaviors for 
the torsional, biaxial proportional and the biaxial non- 
proportional fatigue cycles are consistent with the slip 
deformation analysis discussed in Section 2'. 

4.2.4. Deformation Analyses 

Figures lOa-f show typical optical micrographs of polished 
and etched sections of test samples. These micrographs show the 
orientation of slip bands with refernce to the sample axis (the 
horizontal axis in the micrographs) . Figures 10a and b show the 
micrographs of samples of heat treatment A and F respectively, 
tested in pure torsion. The planes of maximum shear stress in a 
torsional fatigue sample, as shown previously in section 2, are 
parallel and perpendicular to the sample axis. Note that all the 
slip bands are either parallel or perpendicular to the sample 
axis. Figures 10c and d are the micrographs for typical biaxial 
proportional tests. Again, the slip traces seen in these figures 
are oriented approximately at 30° and 120° to the specimen axis, 
and this was as predicted in section 2. Figures lOe and f are 
micrographs for typical non-proportional biaxial tests. As 
predicted in section 2 , the maximum shear planes rotate 
continuously, and therefore, many slip systems become active 
during fatigue cycles. The above descriptions are generally true 
for both heat treatments, but is immediately obvious in heat 
treatment F, because of the larger grain size. 

Figures lla-f show macrophotographs and schematic of sample 
surfaces with the primary cracks. Again, these crack 
orientations are very consistent with the proposed deformation 
modes in Section 2. For example, in pure torsion, the primary 
crack is either perpendicular or parallel to the axis of the 
sample, thus confining to the maximum shear planes (figures 11a 
and b) . Likewise, for the proportional biaxial test, the crack 
is again located on the maximum shear planes (figures 11c and d) . 
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For the non-proportional biaxial tests, the crack is snake shaped 
(figures lie and f) , thus not confining to any specific plane, 
which is consistent with earlier discussions. In all these cases, 
the initial mode II type crack (formed by shear) changes to mode 
I type crack as crack extension takes place. Microcracks on the 
tested sample surfaces were observed in both heat treatments, 
though more frequently in heat treatment A. These microcracks 
were preferentially on the maximum shear planes. Typical example 
of microcracks in a tested sample is shown in figure 12. 
Fractography of the fracture surfaces indicated that the crack 
propagation in these materials is by classical striation 
mechanism. Typical example of the striations is shown in figure 
13. 

5. SUMMARIZING REMARKS AND CONCLUSIONS 

Torsional and biaxial (tension-torsion) fatigue tests of 
Waspaloy were conducted to illustrate the importance of 
metallurgical and microstructural aspects of multiaxial fatigue. 
Waspaloy, in two heat treatments, was tested for its torsional 
and biaxial (tension-torsion) fatigue behavior. The results 
indicate that slip deformation and cyclic hardening are dependent 
upon the type of fatigue test. For example^ biaxial out-of-phase 
(non-proportional) tests were found to be more damaging than 
biaxial in-phase (proportional) tests and torsional tests. This 
is due to the differences in the associated slip deformation. 
Uniaxial, torsional and proportional biaxial fatigue were shown 
to usually exhibit single slip deformation mechanisms and non- 
proportional biaxial fatigue was shown to exhibit multiple slip 
deformation. The slip deformation aspect of multiaxial fatigue 
is important only at low temperatures. At high temperatures, 
other metallurgical effects such as creep damage and 
environmental effects will have significant influence on the 
multiaxial fatigue behaviour. 

The need for defining the test parameters for biaxial tests 
has also been highlighted during the course of this work. For 
example, the conventional definition of plastic strain range, as 
defined by the width of the hysterisis loop for uniaxial fatigue, 
is not any more true for non-proportional multiaxial fatigue. 
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Figure 1. Biaxial fatigue cycles (a) in-phase (proportional) and 
(b) 90° out-of-phase (non-proportional) 

Figure 2. (a) Schematic of a body under biaxial loading, (b) 

Detail of a small section on the surface of the body shown in 

(a) ; accordingly, for (i) pure tension $ = 45° and 135 , (ii) 
pure torsion $ = 0° and 90 , (iii) biaxial proportional $ is 
fixed but dependent on relative magnitudes of A e and A-y , and (iv) 
biaxial 90° out-of-phase 4> continuously rotates in a cycle and 
the rate of change of $ depends on a number of factors. 

Figure 3. Resolved shear stresses for the twelve octohedral slip 
systems in fee crystal during a fatigue cycle, (a) for pure 
tension, (b) for pure torsion, (c) for biaxial proportional cycle 
and (d) for biaxial non-proportional cycle. In figures (a)-(c), 
only one slip system enjoys the privilege of maximum resolved 
shear stress during the complete fatigue cycle, while in (d) many 
slip systems enjoy this privilege, mainly due to the rotating 
maximum shear planes, shown in figure 2. These calculations were 
done for a random orientation. 

Figure 4. Fatigue specimen drawing; all units are in cm. 

Figure 5. High resolution darkfield TEM iamges of the gamma prime 
precipitates, (a) heat treatment A and (b) heat treatment F. 

Figure 6. SEM micrographs of typical grain boundaries, (a) heat 
treatment A and (b) heat treatment F. 

Figure 7. Typical LCF hysterisis loops, (a) torsional fatigue, 

(b) axial loop and (c) torsional loop for in-phase biaxial 
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fatigue, (d) axial loop and (e) torsional loop for 90° out of 
phase biaxial fatigue. 

Figure 8. Cumulative glide plots for (a) torsional fatigue, (b) 
axial stress and (c) torsional stress for in-phase biaxial 
fatigue, (d) axial stress and (e) torsional stress for 90 out- 
of-phase biaxial fatigue. Heat treatment A 

Figure 9. Same as in figure 8, for heat treatment F 

Figure 10. Optical micrographs of polished and etched specimens 
showing slip traces. The specimen axis is horizontal. (a) and 
(b) torsional, (c) and (d) biaxial in-phase, and (e) and (f) 
biaxial 90° out-of-phase for heat treatments A and F 
respectively. Note that the slip traces are parallel to the 
predicted maximum shear planes. 

Figure 11. Macrophotographs showing crack morphology, (a) , (c) 

and (e) are the photographs for torsional, biaxial in-phase and 
biaxial 90° out-of-phase tests and (b) , (d) and (f) are 

corresponding shematic showing the crack orientation. 

Figure 12 . Typical example of optical photograph showing 
microcracks . 

Figure 13. Typical example of SEM photograph showing striations 
on fracture surface* 
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LIST OF SYMBOLS 

r p SS = Resolved shear stress on a slip plane, 

<i> 1 , 4> 2> $ 3 ~ Un it vector in the direction of the slip plane 

normal, 

6 S 2 i 5 3 = Unit vector in the slip direction 
a xx f a yy> a zz = normal stresses at any point 
T xv' T yz> T zx = shear stresses at any point 
a ,r * = applied axial and shear stresses 

$ = angle between specimen axis and plane of maximum 


« / 7 


*i' r o' L 


Ae t, Ae p 

A ?t, A *fp 

N f 


shear 

= applied axial and shear strains 
= angle of twist in a torsional specimen 
= inner and outer radii and gage length of the 
specimen 

= applied torque on the specimen 
= total and plastic axial strain ranges 
= total and plastic torsional strain ranges 
= number of cycles to failure (using the 10%— drop- 
in-load criterion 


13 


Table I 

Composition of 

Waspaloy Plate 

(in Wt. 

Pet. ) 



C 

Co 

Cr 

Mo 

Al 

Ti 

Fe 

Mn 

.046 

13.93 

19.25 

4.16 

1.36 

3.02 

1.17 

. 03 

Si 

S 

P 

Cu 

Sn 

Pb 

B 

Zr 

.05 

.002 

. 006 

. 02 

.0007 

.0001 

. 0051 

.067 


N 

Mg 

Bi 

Ni 

.004 

.0011 

<.250 ppm 

. Balance 


Table II Microstructural Analysis 


Heat treatment A 


Heat treatment F 


Solution Treatment 
Aging Treatment 
Grain Size 


2 hrs. § 1010»C 
Oil Quench 

24 hrs. @ 875 *C 
Oil Quench 

ASTM , #7 

Avg. dia = 321m 


2 hrs. £ 1100*C 
Oil Quench 

6 hrs. @ 730 *C 
Oil Quench 

ASTM #2 

Avg. dia = 180im 


Gamma prime ppte. 
particles size 

Grain boundary 
carbides 


3000 & 900 A 

Numerous, very dense 
distribution 


50 - 100 A 


Coarse grains are free 
of carbides; fine grains 
as in A 


Matrix carbides 


Large blocky carbides Same as in A + fine 

along rolling direction spherical carbides 
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Table III Fatigue Data 


Test 

No. 

Control 

mode 

<f> £ t 
% 

^ « p 
% 

t 

% 

% 

a max a min r max r min N f 

MPa MPa MPa MPa cycles 

F04 

Tor/tot 

_ 

— 

3.00 

1.76 

- 

- 

435 

435 

2,497 


Tor/pla 

- 

- 

2.19 

1.00 

— 

— 

421 

421 

3 , 462 

F07 

Tor/pla 

— 

- 

1.32 

0.30 

— 

— 

372 

372 

10,326 

F09 

Tor/pla 

- 

- 

1.00 

0.10 

— 

— 

138 

138 

21,613 

F01 

InA /pla 

1.74 

1.00 

1.59 

1.00 

717 

772 

228 

228 

1,520 

F02 

InA /pla 

0.91 

0.30 

0.80 

0.30 

648 

676 

- 200 

200 

5,117 

F03 

InA /tot 

0.50 

.088 

1.00 

0.22 

414 

441 

290 

290 

12,573 

Fll 

InA /tot 

0.60 

.060 

0.60 

.074 

572 

600 

207 

214 

11,423 

F08 

OoA/pla 

0.98 

0.30 

1.09 

0.30 

752 

779 

386 

393 

2,439 

F10 

OoA /pla 

1.88 

1.00 

1.98 

1.00 

993 : 

L, 034 

538 

559 

520 

F05 

OoA /tot 

0.50 

.005 

1.00 

0.10 

538 

579 

359 

365 

5,727 

F12 

OoA /tot 

0.60 

. 023 

0.60 

.045 

607 

662 

221 

234 

12,573 

A03 

Tor/tot 

_ 

— 

3 . 00 

1.73 

- 

- 

448 

434 

3,351 

A01 

Tor/pla 

- 

- 

1.38 

0.30 

- 

— 

372 

372 

19,761 

A08 

Tor/pla 

- 

- 

2.41 

1.00 

- 

— 

434 

428 

3,705 

A09 

InA /pla 

1.03 

0.30 

0.87 

0.30 

676 

724 

207 

241 

7,311 

All 

InA /pla 

1.92 

1.00 

1.69 

1.00 

731 

841 

228 

221 

1,429 

A10 

InA /tot 

0.50 

.064 

1.00 

0.18 

434 

483 

310 

296 

24,214 

A12 

InA /tot 

0.60 

.054 

0.60 

. 063 

545 

607 

207 

207 

33,350 

AO 2 

OoA /pla 

1.10 

0.30 

1.20 

0.30 

786 

821 

386 

400 

2,620 

AO 6 

OoA /pla 

2.00 

1.00 

1.69 

1.00 1,014 : 

L , 048 

503 

524 

429 

A05 

OoA /tot 

0.50 

.007 

1.00 

.090 

503 

579 

345 

352 

11,029 

AO 7 

OoA /tot 

0.60 

.026 

0.60 

.046 

600 

669 

221 

228 

24,218 


Notes: 

(i) Control modes: Tor = Torsion; InA = in-phase; OoA = 90* out-of-phase 

tot = total strain range controlled 
pla = plastic strain range controlled 

(ii) ^ t = total axial strain range; 

<f> e p = axial plastic strain range (axial fatigue loop width); 

= total tosional strain range; 

<f, $ p = torsional plastic strain range (torsional fatigue loop width) 
a max = max i muin axial stress (tensile) ; 

^min = m i n i- muin axial stress (compressive) ; 
r = maximum torsional stress; 

luaX 

r min = torsional stress; 

= number of cycles to failure. 

(iii) All of the above parameters are sampled at N^/2. 

(iv) Fatigue loop widths do not really indicate the plastic strain range 
for the 90* out-of-phase biaxial fatigue tests, since for these tests 
the elastic unloading during a cycle is incomplete or totally absent; 
those data are however included here for the sake of completeness. 
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Torsional Strain 


(a) Axial-torsional LCF (in- 
phase) axial loop 



(c) Axial-torsional LCF (in 
phase) torsional loop 
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loop 


(e) Axial-torsional LCF (90' 
out-of-phase) torsional 
loop 
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Figure 7 


Torsional Strain 



Axial or Shear Stress 


(a) Torsional LCF ; plot of shear stress 
vs. cycles. 


(b) Axial-torsional (in-phase) LCF 
plot of axial stress vs. cycle 


(c) Axial-torsional (in-phase) LCF; 
plot of shear stress vs. cycles 


(d) Axial-torsional (90° out-of-phase) 
LCF; plot of axial stress vs. cycles. 


(e) Axial-torsional (90° out-of-phase) 
LCF; plot of shear stress vs. cycles. 


Number of Cycles 

Figure 8 







Axial or Shear Stress, 



(c) Axial-torsional (in-phase) LCF; plot 
of shear stress vs. cycles. 



(d) Axial-torsional (90° out-of-phase) 
LCF; plot of axial stress vs. cycles 


(e) Axial-torsional (90° out-of-phase) 
LCF; plot of shear stress vs. cycles 


Number of Cycles 


Figure 9 
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BLANK NOT FILMED 

A life prediction method for tension-torsion fatigue: 

ABSTRACT 

A life prediction model has been developed using a cyclic 
equivalent strain approach for predicting biaxial (tension - 
torsion) fatigue lives based on uniaxial fatigue data. The 
cyclic equivalent strain approach has enabled us to combine 
uniaxial, torsional and biaxial (tension-torsion) fatigue data 
into a single life equation. The above model was applied 
successfully to experimental data available in Waspaloy (two heat 
treatments) and 1045 steel. A comparison with other life models 
has also been made. 

1 . INTRODUCTION 

Many strucural members and machine components are subjected 
to repeated loading and therefore fatigue is one of the major 
considerations in their design. A number of life prediction 
methods are available for uniaxial cyclic loading (1). Actual 
structural elements, however, may often be subjected to 
multiaxial cyclic loading and a life prediction methodology 
common to these different loading conditions is essential for 
design applications. Once such a method becomes available, one 
can predict multiaxial fatigue life from uniaxial fatigue data. 
Most of the current multiaxial fatigue life prediction criteria 
(2) are limited to proportional loading conditions* and 
application of these criteria to non-proportional cyclic loading 
leads to significant error. 

This paper describes a Cof f in-Manson type fatigue life 
prediction methodology for complex biaxial loading conditions 
based on cyclic equivalent strain. Under proportional loading 
conditions, the cyclic equivalent strain is determined by using 
Von-Mises criterion and under non-proportional loading 
conditions a correction factor to the equivalent strain is 
introduced based on the strain path, which is phase angle 
dependent. Through this correction factor it is shown that one 
can normalize any type of loading into cyclic equivalent strain 
and the method is shown to apply effectively to experimental data 
available for Waspaloy and 1045 steel. 

2. EQUIVALENT STRAIN APPROACH 

Equivalent strain has been previously defined (3) for in- 
phase biaxial fatigue as: 
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♦Note : Proportional loading refers to fatigue cycles in which 
the axial and torsional cycles are in-phase; non-proportional 
refers to conditions when they are not in-phase. The differnces 
are illustrated in figure 1. 
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2 ,h 


e = 


[ ( £ i~ e 2 ' + ( € 2~ e 3 ' + (e 3 -e l) ^ 


( 1 ) 


where e e o an< ^ f 3 are principal strains. 
When condition of constancy of volume, i.e., 

€ l + e 2 + € 3 = ® * 

is invoked, equation (1) reduces to 


= {- I - [(« i ) 2 + (< 2 ) 2 + 


(e 3 ) 2 ])^ 


( 2 ) 


"3 — l v" p ~ 2 ' '’3- 

When a thin walled tubular specimen is subjected to axial 
strain, e and torsional strain, c simultaneously, the applied 
strains are represented by the following strain tensor. 


[7] = hi -he 0 

0 0 -he 

The equivalent strain e is obtained by calculating the 
principal strains of e , and then substituting these into 
equation (2) 


r = {e 2 + ( 7 2 / 3 ) } h * • • (4) 

Strain states leading to similar e values, regardless of the 
individual values torsional and axial strains, are expected to 
result in similar fatigue lives. Use of the above equivalent 
strain is justified only when the directions of the principal 
strains do not change during cyclic loading ( 4 ), i.e., under 
proportional loading conditions. For non-proportional cyclic 
loading, the principal axes rotate continuously and elastic 
unloading is normally incomplete. 

The net effect of rotation of the principal axes and the 
elastic unloading behaviour for non-proportional cyclic loading 
is the increase in the strain path per cycle and the increase in 
the area under the e - 7 strain space. For example, as shown in 
figures 2a and b, the strain path per cycle for non-proportional 
loading is ! 2 n(e a 2 + 7 2 )' i , while for proportional loading it is 
4 ( c 2 + 7 2 )^. Thus, the strain path for non-proportional 
loafing is 11.1% greater than that for proportional loading. 

The damage due to this increased strain path will depend on 
the amount of elastic unloading (or the lack of it) , which in 
turn depends on the phase angle, <j> , between the axial and the 

torsional strain. If one assumes that the area under the e - 7 
space is indicative of the elastic unloading, i.e., the area is 
inversely proportional to elastic unloading, a correction factor, 
f(^), for the cyclic equivalent strain given in equation ( 4 ) can 
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be developed. It is easily shown that the area under the e - 7 
space is proportional to Sin <f> . 

Combining the increase in the strain path per cycle and the 
area under the e - 7 space, the correction factor becomes: 

f(^) « 0.111 Sin <j> t ••• (5) 

Now the generalized cyclic equivalent strain is given by: 

€ cycle = U + f (*>> t € a + ^a ^ ••• ( 6 ) 

where the correction factor takes a value of zero for 

proportional loading and a maximum value of 0.111 for 90 out-of- 
phase, non-proportional loading. For intermediate phase angles, 
the correction factor falls between these two extreme values. 

The underlying idea of this simple concept is that for a 
given set of axial and torsional strain ranges in a biaxial test, 
an increase in non-proportionality will result in considerably 
increased hardening and therefore it has the same effect as 
increasing the applied strains for a proportional biaxial test. 
For example, when in-phase and 90° out-of phase biaxial fatigue 
tests were conducted under similar axial and torsional strain 
ranges, the out-of-phase loading was found to result in much 
higher stabilized stresses and shorter fatigue lives (5,6). 

The generalized cyclic equivalent strain given in equation 
( 6 ) can then be used in a Cof f in-Manson type life relationship of 
the form: 

Nf = A (e C y C ie) *'* ( 3 * * * 7 ^ 


3. VERIFICATION OF THE MODEL 

The above life prediction method enables us to combine 

experimental data from uniaxial fatigue tests, torsional fatigue 
tests and proportional and non-proportional biaxial fatigue 
tests. In this section, the model is applied to experimental 
data from Waspaloy (two heat treatments) and 1045 steel. 

3.1. Waspaloy 

Lerch, Jayaraman and Antolovich (7) had conducted uniaxial 
fatigue tests on Waspaloy. In their study, Waspaloy was 
subjected to two heat treatments, referred to as A and F in 
reference ( 7 ) , resulting in microstructures of small grain size 
(ASTM grain size 9) and coarse gamma prime precipitate particles 
(about 900 A in diameter) for heat treatment A and large grain 
size (ASTM grain size 3) and fine gamma prime precipitate 
particles (about 50-80 A in diameter) for heat treatment F. In a 
follow-up study, torsional and biaxial tension-torsion (in-phase 
and 90° out-of-phase) fatigue tests for similar microstructural 
conditions were conducted by Ditmars and Jayaraman (5). The 
experimental results from these two studies were analyzed using 
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the cyclic equivalent strain approach given in equations (6) and 
(7), and the corresponding Cof f in-Manson type plots are shown in 
figures 3(a) and (b) . As seen in these two figures, the data 
from uniaxial and multiaxial fatigue tests can be represented in 
one life curve using the cyclic equivalent strain approach. 
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3.2. 1045 Steel 

Fatigue test data were obtained from two sources. Appel and 
Jayaraman (6) conducted in-phase and 90° out-of-phase tension- 
torsion tests on 1045 steel. Fash et al. (8) conducted uniaxial, 
torsional and biaxial in-phase fatigue tests on 1045 steel. Data 
from these two sources were analyzed using the cyclic equivalent 
strain approach and the resulting life plot is shown in figure 4. 
Again this figure shows that data from uniaxial and multiaxial 
fatigue tests can be represented on one life curve. 

3.3. Life prediction 

In another exercise, attempts were made to predict biaxial 
in-phase and out-of-phase fatigue lives using only the uniaxial 
fatigue data for Waspaloy and 1045 steel. Here the uniaxial 
data for Waspaloy from reference 7 and that for 1045 steel . from 
reference 8 were used to get the constants A and b in equation 
(7) . Using these constants, fatigue lives for biaxial conditions 
(both proportional and non-proportional) corresponding to the 
experimental data for Waspaloy from reference 5 and that for 1045 
steel from references 6 and 8 were predicted. Figure 5 shows the 
plot of predicted life vs. actual life. As- seen in this figure 
most of the predicted data fall within the factor of two limits. 


4. COMPARISON WITH OTHER LIFE PREDICTION METHODS 

Zamrick and Frishmuth (9) used a "total strain", e T , 
defined as: 


/ 2 I 2 I 2\h 

e T -( e l + € 2 + € 3 ^ 


... ( 8 ) 


for correlating life. For biaxial case e T is calculated using 
the equation given below: 


« ip = { [ \ ([e" x «t — v e " x Sin wt] 


+ { ( e " Y Sin wt + i/e» y Sin wt ) 2 + (7 " xy Sin (wt+^ ) ) » } *) ] 2 


x w w ' x 
+ [ ^ ([« " v Sin wt - v e " Sin wt] 


2 \ h \ i 2 


- {(e"x sin wt + 1/6 » x Sin wt)^ + (7 " xy Sin (wt+^ ) ) *) *») ] 


+ [-»/e " x Sin wt]'* } 


2 x h 


... (9) 


where , 

« " x » 7 " xy = maximum applied axial and torsional strains. 


w 


= frequency of straining, 
= Poisson's ratio and 
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t = time 

(emJjnax was computed numerically using the above equation for 
each set of experimental data. Figures 6 a and b plot (« vs. 
log Nf for the fatigue data in Waspaloy in the two heat 
treatments A and F respectively. Figure 6 c shows a similar plot 
for the fatigue data in 1045 steel. Data correlation in these 
three figures are good. If the non-proportional data are 
disregarded, most of the other data points lie within a factor of 
2 on life. The theory is unfortunately non-conservative for out- 
of— phase tests. Also, one of the major drawbacks of this model 
is that the calculation of («m) m -x involves determining the 
maximum of the function described in equation 9 numerically and 
these calculations are even more complex for phase angles other 
than 0 ° or 90°. 


Garud (10) used inelastic work (W ) per cycle as a measure 
of damage. For a biaxial test, inelastic work is simply the sum 
of axial and torsional hysterisis loop areas. Figures 7a and b 
plot log W D vs. log N f for fatigue data in Waspaloy in its two 
heat treatments A and F respectively. Again as seen in these 
figures, the data correlation is not good since many of the data 
points are outside the factor of 2 on life. In addition, one 
major problem in using this method is that prediction of Wp 
without actually running an experiment is very difficult. 


The method proposed by Socie and Shield (11) is 
representative of a class of theories that are concerned with 
forces acting on some specific plane, rather than global 
quantities, which the previous methods considered. The. main 
problem with these theories is that they can only be applied to 
proportional data. Therefore, non— proportional data will not be 
considered in the evaluation of this method. Socie and Shield 
use the quantity 


A A A 

D = 7 p + 1.5 e np + 1.5 r n /E ••• (10) 

as a measure of damage. Figures 8 a and b plot D vs. log Nf for 
fatigue data from Waspaloy in the two heat treatments A and F 
respectively. As indicated earlier, the non-proportional data 
were not included and inspite of this, as seen in these figures, 
the data correlation is not very good. 

5. SUMMARY 

A fatigue life model using cyclic equivalent strain approach 
has been developed to predict fatigue lives under multiaxial 
fatigue conditions. This model has been shown to apply well to 
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experimental data from Waspaloy and 1045 steel. As a comparison, 
three earlier models were applied to the same data base and the 
deficiencies in unifying all fatigue data have been indicated. 
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Figure 1. Biaxial fatigue cycles (a) in-phase (proportional) and 
(b) 90° out-of-phase (non-proportional) . 

Figure 2. e - j space for strain controlled multiaxial fatigue, 
(a) in-phase and (b) 90° out-of-phase. 
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Figure 3. Unified fatigue life plots for Waspaloy, (a) heat 
treatment A and (b) heat treatment F. 

Figure 4. Unified fatigue life plot for 1045 steel. 

Figure 5. Plot of predicted life vs. actual life for torsional 
and biaxial fatigue; the predictions are based on constants 
obtained from uniaxial fatigue data using the cyclic equivalent 
strain approach. Data from two heat treatments of Waspaloy and 
1045 steel are presented. 

Figure 6. Fatigue life plots using Zamrick and Frishmuth method, 
(a) Waspaloy, heat treatment A, (b) Waspaloy, heat treatment F 
and (c) 1045 steel. 

Figure 7. Fatigue life plots for Waspaloy using Garud method, (a) 
heat treatment A and (b) heat treatment F 

Figure 8. Fatigue life plots for Waspaloy using Socie and Shield 
method, (a) heat treatment A and (b) heat treatment F 

NOMENCLATURE 


e 1 , e 2 ' e 3 = P r i c iP a l strains 
e ,7 = normal and shear strain 


£ a /7 a = axial and torsional strain amplitudes 

£ = von Mises' type equivalent strain when axial and 

torsional cycles are in-phase 


€ cycle 


cyclic equivalent strain under all biaxial conditions 


= phase angle between axial and torsional cycles 


f (*) 


N f 

A and b 


€ 


ii ~ it 

x' 7 xy 


e T 

CO 

t 


u 

D 



7 P 
a n 


= correction factor to e for phase angles v = 0°; f(v) 
takes values between 0 for v = 0° and 0.111 for v = 90° 
= cycles to failure 
= constants in life equation 

= maximum applied axial and torsional strains as defined 
in reference 9 

= total strain as defined in reference 9 
= frequency of straining 
= time 

= Poisson's ratio 
= plastic damage parameter 

= plastic axial strain amplitudes normal to maximum 
shear planes (reference 11) 

= plastic shear strain amplitude on the maximum shear 
plane (reference 11) 

= stress amplitude normal to maximum shear plane 
(reference 11) 
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Nf = 0.002 ( £ Cycle ) 


R = -0.982 


factor of two 


* margin 
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